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The two-dimensional spin-gap system SrCu2(B03)2 shows unique physical properties due to the 
low-dimensionality character and the strong quantum fluctuations. Experimentally, 1/8-, 1/4-, and 
1/3- plateaus have been observed in the magnetization curve under magnetic fields up to 70 Tesla, 
and the 1/2-plateau is expected to be stabilized at higher magnetic fields. We argue that spin- 
lattice effects are necessary to describe the superstructures at the plateaus, and we propose a simple 
microscopic model of spins interacting adiabatically with the lattice to reproduce the main features 
of the recent experimental results by nuclear magnetic resonance. 

PACS numbers: 75.10.Jm, 75.30.Kz 



I. INTRODUCTION 

The possibility to obtain low-dimensional quantum 
spin systems that do not order magnetically at very low 
temperature, or even down to zero temperature, is cur- 
rently a subject of great interest. Although there are 
many examples of onc-dimensional (ID) or quasi-lD sys- 
tems such as SrCu 2 3 (5=1/2 ladder),! Y 2 BaNi0 5 
(Haldane chain) j 2 . CuGeC>3 and L1V2O5 (prototype of 
frustrated S — 1/2 chains),^ for some time, the only 
example of a two-dimensional (2D) system with a singlet 
ground state and a finite gap to magnetic excitations 
was the vanadium oxide CaV40g (Ref. 000)- There- 
fore, the recent discovery of SrCu2 (603)2 (Ref. @) rep- 
resents a breakthrough in this direction and paves the 
way for further research on this unconventional state of 
matter. This compound has a layered structure where 
stacking layers of CUBO3 are intercalated by magnet- 
ically inert layers of Sr. A spin S = 1/2 resides on 
each Cu 2+ ion, forming a 2D orthogonal dimer lattice^ 
(see Fig. ^| . It turns out that the magnetic properties 
of SrCu2(BC>3)2 are very well described by the 2D or- 
thogonal dimer Hciscnberg model^ which is topologically 
equivalent to the 2D Shastry- Sutherland models 
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FIG. 1: The orthogonal dimer lattice. The continuous and 
dashed lines indicate the antiferromagnetic interactions and 
the elastic couplings for nearest neighbor and next-nearest 
neighbor sites, respectively. The nearest neighbor and the 
next-nearest neighbor equilibrium distances are indicated by 
do and d' , respectively. The 16-site square cluster and 24-site 
rectangular cluster used in the calculations are shown by thin 
solid lines. The transferred hyperfine couplings B and C (see 
Section HV1 are also shown. The arrows define the direction 
from the site i to the site j for the Dzyaloshinsky-Moriya 
interaction (see Section flV Bl . 



here Si = (Sf, Sf , Sf) is the spin-1/2 operator at the site 
i and the notations (n.n.) and (n.n.n.) stand for nearest 
neighbor and next-nearest neighbor sites, respectively. In 
the parameter range J'/ J < 0.68 (Ref. Hi"!) the ground 
state is exactly known to be the product of dimer sin- 
glets J£ 
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where a indicates the dimer bond connected by the su- 
perexchange coupling J. In Ref. ^2 it has been shown 
that the values of the antiferromagnetic superexchangc 
are such that J'/ J = 0.635, with J — 85 Kelvin, indi- 



cating that SrCu2(B03) 2 can be described by a 2D spin 
system whose ground state is exactly known. 

The very particular nature of the ground state induces 
unique features in the spin excitations at low temper- 
atures. First of all, there is a finite gap to the mag- 
netic excitations. This spin gap has been observed in 
several experiments and was estimated to be about 35 
j£ e j v j n( i3 i i4 i i5 i i6 ]vi ore0 ver, the triplet excitations have 
an almost localized nature because of the orthogonality 
of the J bonds^ Such a behavior was revealed experimen- 
tally, using inelastic neutron scattering^ as an almost 
flat triplet dispersion. 

The localized nature of the triplet excitation leads 
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also to the most spectacular phenomenon of this sys- 
tem: When an external magnetic field, up to 70 Tesla, is 
applied, the system shows magnetization plateaus corre- 
sponding to 1/8, 1/4, and 1/3 of the full Cu 2+ moment 
(another plateau at 1/2 of the full Cu 2+ moment is likely 
to be stabilized at even higher mag netic fields) ML By 
applying an external magnetic field, the density of the 
triplets can be tuned and it is found that the magneti- 
zation stays constant for particular ranges of the exter- 
nal field. It is expected that the plateaus are due to a 
crystallization of the triplets for particular commensu- 
rate values of the magnetization and are due to the small 
ratio between the kinetic and the interaction energy of 
the triplets. To our knowledge, SrCu2(B03)2 is the first 
example of a 2D quantum spin system which shows mag- 
netization plateaus. 

So far, the magnetization plateaus have been studied 
mostly in ID systems ^i^iSSiSi In particular, Yamanaka, 
Oshikawa, and Affleck 2 ^ have found a simple necessary 
condition for the existence of plateaus in ID systems. De- 
noting by I the period of the ground state in the presence 
of the external field, by 5 the magnitude of the spin and 
by m the magnetization per site (in unit of .g/is); the oc- 
currence of a magnetization plateau is only possible when 
the condition 

l(S — m) = integer, (3) 

is satisfied. It is important to mention that I can be 
different from the period of the lattice. The first exam- 
ple of a magnetization plateau accompanied by such a 
symmetry breaking is a 5 = 1/2 Heisenberg chain with 
next-nearest-neighbor and alternating nearest-neighbor 
interactions 

Recently, Oshikawa extended the Licb-Schultz-Mattis 
argument to quantum many-particle systems with a con- 
served particle number on a periodic lattice in arbitrary 
dimensions and showed that the condition of Eq. for 
the plateau is still valid in arbitrary dimensions^ In- 
deed, several theoretical works concerning the plateaus 
for SrCu2(BOa)2 indicate that the criterion of Eq. 10 is 
satisfied and, except for the plateau at 1/2, the ground 
states at the plateaus are accompanied by a breaking of 
the translational symmetryi 26 i 27 i 28 i 29 i 30 i 31 

Recently, Kodama and collaborators 3 ^ performed a nu- 
clear magnetic resonance (NMR) measurement at the 
1/8-plateau, for an external magnetic field H = 27. 6 
Tesla and a temperature of 35 mK, and observed the 
presence of at least 11 different sites (i.e., different values 
of the local magnetization), indicating a clear breaking 
of the translational symmetry. So far, from a theoretical 
point of view, the superstructures at the plateaus have 
been studied only by an effective hard-core boson model. 
In this approximation, the triplet with S z = 1 is repre- 
sented by a hard-core boson and the dimer singlet by a 
vacancy. An effective Hamiltonian is derived by pertur- 
bation theory and the magnetization curve and super- 
structures have been calculated by solving it, 2 iii2L22i22i21 
In this picture, two different unit cells with 16 sites have 



been proposed to describe this plateau: i) a 16-site square 
unit cell and ii) a rhomboid unit celli^ However, the 
ground state at the 1/8-plateau is described by a state 
where one of the eight singlets is promoted to a triplet 
within the unit cell. Thus there are only two different 
sites, corresponding to the singlet and the triplet states, 
and it cannot reproduce the very rich texture of the mag- 
netization observed in NMR experiments. Therefore, in 
order to reproduce the experimental data, the simple 
hard-core boson model is not sufficient and it is neces- 
sary to consider the original spin Hamiltonian. 

From the results of the hard-core boson calculations, as 
well as from the result of Ref. [2f|, it is expected that the 
ground states at plateaus are degenerate. On the other 
hand, a magnetization structure will only be observed 
if an extra mechanism selects one of the ground states, 
since linear combinations of the ground states can lead 
to an arbitrary magnetization. For instance, a uniform 
linear combination would lead to a uniform magnetiza- 
tion. In the real material, this mechanism could be due 
to pinning by impurities, or to a lattice distortion. Ac- 
tually, a strong motivation for considering lattice effects 
in SrCu2(B03)2 is given by the pronounced softening of 
the sound velocity observed at the edges of the magne- 
tization plateaui 3 ^*^ On the theoretical side, one could 
try to select a magnetization texture by imposing an ex- 
ternal, very small symmetry breaking field, as is often 
done for instance for the dimerized state in a spin-Peierls 
system. However, in the present case, we do not know 
a priori the magnetization texture, and imposing a spe- 
cific field would bias the results. Therefore, we prefer the 
more physical way which consists in coupling the system 
to phonons. 

The paper is organized as follows. In Section [HJ we 
introduce the spin-phonon Hamiltonian and explain the 
method. In Section InTl we present the results for the su- 
perstructures at 1/8-, 1/4-, 1/3- and 1/2-plateaus. In 
Section IIVI we compare our theoretical results for the 
1/8-plateau to the experimental ones, including in addi- 
tion the effects of inter-layer coupling and Dzyaloshinsky- 
Moriya interaction. Finally, Section [V] is devoted to the 
conclusions and the discussion. 



II. MODEL AND METHOD 

In this Section we introduce the spin-phonon Hamil- 
tonian and we describe the method that we use to 
characterize the spin texture of the different magneti- 
zation plateaus. We consider the 5=1/2 orthogonal 
dimer model coupled to adiabatic phonons that we study 
by exact diagonalizations of finite clusters with a self- 
consistent Lanczos algorithm. In this approach, the adi- 
abatic phonons are described by classical variables, re- 
lated to the displacements of the lattice sites. The full 
spin-phonon Hamiltonian is defined on a 2D orthogonal 
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dimcr lattice of N sites by: 



£ ^'(<i„)s.-s. ■ K '(V*- S *> 

(n.n.n.) 




(4) 



here J{dij) and J'(dij) are the antifcrromagnctic su- 
perexchange couplings, which depend on the relative dis- 
tance dij = ||R°+5rj — R° — £rj|| between sites i and j. K 
and if' are the elastic coupling constants, Sdij = dij—cqj, 
and d?- = |jR°-R°|| is the cquilibr ium distances between 
Copper sites. 

For small displacements of the Cu sites, it is possible 
to linearize the antiferromagnetic couplings around their 
equilibrium values, Sdij ~ (R" — R°) • (Sri — Stj), and 
therefore we expect that in general: 35 
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In the following, we will denote by a and a 1 the two 
(in principle different) exponents for J(dij) and J'(dy), 
respectively. 

With a self-consistent Lanczos diagonalization of finite 
clusters, one can find the optimal configuration of the 
bond lengths and local spin configurations for given val- 
ues of the coupling parameters (a, a', K and K') and 
given total spin S quite easily: Starting from a random 
choice of the atomic displacements, one just has to im- 
prove iteratively the total energy by changing the lattice 
parameters until a stationary configuration is reached. In 
all cases, we have verified that this configuration is sta- 
ble and unique (up to trivial symmetry operations) by 
starting from different initial distributions, so it must be 
the global minimum. This method has been already used 
for other ID, quasi-lD and 2D spin systems interacting 
adiabatically with the lattice^SiSi With respect to other 
approximate approaches, this method has the great ad- 
vantage that it gives unbiased results even for strongly 
frustrated systems, where other numerical methods can 
be highly questionable. It is worth noting that, in our 
self-consistent Lanczos method, the only approximation 
is to consider the adiabatic limit for the phonons. Once 
we restrict our calculation to this case, we obtain ex- 
act results for the lattice distortions on the chosen finite 
cluster. 

In the following, we consider unit cells containing 
N = 16 sites or, in some case, up to N = 24 sites with 
periodic boundary conditions (see Fig. ^| . In our simple 
microscopic model, we assume that the magnetic field 
directly couples to the total spin of the system, stabi- 
lizing the states with higher total spin. Because of the 
finiteness of our cluster, in order to study the plateau at 
1/n, we assume that the external field is able to stabilize 
a state with a given magnetization and we fix the total 



spin of the cluster (by fixing S% ot = J2i = It is 
worth mentioning that, within our self-consistent Lanc- 
zos method, we are not able to study the actual width of 
the magnetization plateau in the thermodynamic limit. 
Nonetheless, once the existence of a given plateau is as- 
sumed, we can produce important insight into the local 
structure of the spins inside the cell. 

In the following calculations, we consider the case of 
a' = 1.75a and K = K' = 750J' but we want to stress 
that all the results are not specific of this particular 
choice and, qualitatively, we found similar results also 
for different values of these parameters. The parameters 
K and K' are only effective parameters that cannot be 
directly matched to the phonon dispersion. An appropri- 
ate model of elastic constants should include the springs 
between all nearest neighbor atoms, but this does not 
improve the calculation because the dependence of the 
spring constants on the actual position of all atoms is 
not exactly known. Nevertheless, the order of magnitude 
of the elastic constants is expected to be similar in all 
oxides. For instance, for CuGc03, in order to have rela- 
tive displacements smaller than a percent - as observed 
experimentally in the spin-Pcicrls phase - the values of 
the elastic constants are required to be of the order of 
100000 Kelvin^ Finally, superexchange theory suggests 
that typical value for a in oxides are in the range 6-14 
(Ref.E3). 

The bond length of nearest neighbor bond do is taken 
as the unit of the length, which is 2.9lA at 100 Kelvin, 
and the length of next-nearest neighbor bond d' is as- 
sumed to be given by d' Q = 1.75d (Ref. l39|) . Finally, 
we note that the Hamiltonian 10} is invariant under the 
rescaling a — > Xa, a' — > Xa', K — ► A 2 it, K 1 — ► A 2 if', and 
Sri ^ <fa*j/A, where A is the rescaling parameter. This 
allows us to fix one parameter among a, a', K and K' . 
Note that the physical values of the magnetic couplings 
are unaffected by this transformation. 



III. EFFECT OF SPIN-PHONON COUPLING 
AT PLATEAUS 

A. 1/8-plateau 

In order to describe the two most probable triplet pat- 
terns that were previously suggested for the 1/8-plateau 
by the hard-core boson approach^ we have diagonalized 
the Hamiltonian Q on two 16-site clusters correspond- 
ing to different unit cells, the 16-site square cluster and 



the rhomboid cluster in the sector with Si, 
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Figs. |2 and The results for the two cases turn out to 
be qualitatively different: In the square cluster we find 
that the ground state has only six different sites, that is 
six different values of the local magnetization (S* ) (see 
Fig. I2J , whereas the rhomboid cluster contains eight dif- 
ferent sites with different local magnetization (see Fig.|3J). 
The energy difference between these two states is very 
small (of the order of 10~ 5 -j- 10~ 6 J per site). In both 
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FIG. 2: Spin density profile for 1/8-plateau for the 16-site 
square cluster. Full (empty) circles indicate sites with mag- 
netization along (opposite to) the external field and the size 
of the circles is proportional to the spin amplitude. The bond 
lengths for a = 10 are also shown. 



cases, the magnetization is centered around one strongly 
polarized dimcr, that, in the following, we will denote 
by "triplet" , with Friedel-like oscillations in the spin am- 
plitude. It is worth noting that we obtain both positive 
and negative magnetizations, and in particular, in both 
clusters, there is a large negative spin, just near to the 
strongly polarized dimer. 

One of the main features of these results is the exis- 
tence of two sites with large and positive polarization and 
one site with large and negative polarization, which is in 
agreement with the NMR experimental finding^ The 
Friedel-like oscillations decay quite fast in space, and far 
from the "triplet" dimer the local magnetization is very 
small. For this reason, we expect rather small size effects 
for our finite cluster calculation and we believe that even 
the 16-site lattice can represent quite well the real sys- 
tem. Of course, we cannot rule out that a larger number 
of different sites (with very small magnetizations) exists 
when a bigger cluster is considered, but the very fast de- 
cay in the oscillation of the local magnetization clearly 
indicates that our small 16-site cluster is able to capture 
the main ingredients of the true ground state: Two sites 
with large and positive (Sf ), one site with large and nega- 
tive (Sf ) and a bunch of sites with rather small (positive 
and negative) magnetizations. Unfortunately, the fact 
that we do not want to impose any external constraint 
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FIG. 3: Spin density profile for 1/8-plateau for the rhomboid 
cluster of 16 sites. Full (empty) circles indicate sites with 
magnetization along (opposite to) the external field and the 
size of the circles is proportional to the spin amplitude. The 
bond lengths for a = 10 are also shown. 

on the spin pattern prevents us from using the lattice 
symmetries in the Lanczos diagonalization, and, in order 
to obtain the optimal lattice (and spin) configuration, it 
is necessary to perform the diagonalization several times 
up to convergence. These two facts make the calculation 
very heavy, and the next cluster, with 32 sites, which is 
consistent with the 1/8-plateau, is beyond the present 
computational possibilities. 

In the optimized lattice configuration, all the transla- 
tional symmetries are broken. The ground state is eight- 
fold degenerate for the square unit cell. On the other 
hand, in the rhomboid case, two types of unit cell arc pos- 
sible and so the ground state is sixteen-fold degenerate. 
The behavior of the local magnetization as a function of 
the spin-phonon coupling is reported in Figs. 0] and [S] for 
the square and the rhomboid cluster, respectively. Notice 
that below a critical value of the spin-phonon coupling a, 
a state with a different number of sites is stabilized (four 
sites for the square cluster and only one for the rhom- 
boid cluster). This is due to the fact that in the 16-site 
square (rhomboid) cluster without spin-phonon coupling, 
the ground state is four-fold (two-fold) degenerate, and 
the other low-lying states are separated by a finite-size 
gap. Thus for small enough spin-phonon coupling, the 
ground state corresponds to a linear combination of only 
four (two) states, and it is only possible to mix all the 
eight low-lying states with a sufficiently large coupling, 
resulting into a state with six (eight) different sites. The 



16-site square cluster for the 1/8-plateau as function of the 
spin-phonon coupling a. 
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FIG. 5: Local magnetization for the eight different sites of 
the rhomboid cluster for the 1/8-plateau as a function of the 
spin-phonon coupling a. 
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FIG. 6: Spin density profile for 1/4-plateau. Full (empty) 
circles indicate sites with magnetization along (opposite to) 
the external field and the size of the circles is proportional 
to the spin amplitude. The bond lengths for a = 10 are also 
shown. 



lengths for a = 10 for square (rhomboid) unit cell are 
shown in Fig. [21 (Fig. |2J. However, for realistic values of 
the elastic constants and the spin-phonon coupling, these 
displacements are very small - of the order of one percent 
or less - and we do not think that there is any chance 
at the moment to detect such a tiny structural distortion 
with X-Rays in a field of 27.6 Tcsla. 



B. 1/4-plateau 



critical value of the spin-phonon coupling is much larger 
for the rhomboid cell (of the order of 5) than for the 
square unit cell (between 1 and 2), see Figs. 01 and |SJ As 
in the spin-Peierls case^ it is expected to decrease by 
enlarging the lattice size; unfortunately, for reasons given 
just above, it is impossible to attempt a size scaling of 
this value. 

For both unit cells, we recover the hard-core boson 
results of Ref. [2^ in the extreme limit of infinite spin- 
phonon coupling (a — » oo), where all the magnetization 
is carried by a localized triplet dimer, and all the other 
dimers are perfect singlets. 

Finally, we want to make a remark on the lattice dis- 
placements. The magnetic energy gain is related to a 
shrinking of the bond length of the "singlets", whereas 
the "triplet" dimers enlarge their bond length. The bond 



Experimentally the next magnetization plateau is at 
1/4 of the total Cu 2+ moment. From the hard-core boson 
calculations^ it comes out that the unit cell is just half 
of the previous 16-site square cluster. Therefore, for this 
plateau, we consider the 16-site square lattice and we 
perform the self-consistent Lanczos method in the sector 
Sf ot = 2. Moreover, we also report some results for a 
larger 24-site rectangular lattice with S^ ot = 3. 

The typical outcome of our exact calculation is shown 
in Fig.El Wc find four different sites, three with a positive 
and one with a negative magnetization. Also in this case, 
the "triplet" dimer is almost localized and the system 
shows Friedel-like oscillations on the spin amplitude, the 
nearest site to the "triplet" having negative (S z ). The 
ground state is eight-fold degenerate and the eight states 
are found to be connected by simple symmetry operations 
of the lattice, like translations and/or reflections. The 



the 1/4-plateau as a function of the spin-phonon coupling a. 
The full dots indicate the results for the 16-site cluster and 
the empty squares the results for the 24-site one at a = 10 
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behavior of the local magnetization as a function of the 
spin-phonon coupling is reported in Fig. \7\ and it clearly 
indicates that by increasing a, the total magnetization 
concentrates progressively in the "triplet" dimer and the 
local magnetizations of the other dimers tend to zero. In 
the extreme limit of a — > oo, the picture corresponds to 
the hard-core boson approximation, where there is only 
one triplet that carries all the magnetizationi^SiSl 

Notice that the spin texture is built up by diagonal 
stripes of strongly polarized dimers, i.e., "triplets", in- 
tercalated by dimers with a small polarization. This pic- 
ture comes out from the fact that there are repulsive 
interactions between the triplets, which make the stripe 
structure stable (this can be seen for example in the per- 
turbation theory 28 ! 30 ). 



C. 1/3-plateau 

In order to study the 1/3-plateau, we are forced to con- 
sider the 24-site cluster with S* ot = 4 because in the 16- 
site cluster the corresponding spin sector is not present. 
It is worth noting that this cluster breaks the reflection 
symmetries of the original orthogonal dimer lattice; on 
the other hand, the shape of this cluster fits well the 
pattern suggested in previous hard-core boson calcula- 
|.j onSi 26 i 27 i 28 J 29 p or .j- n j s pj a teau, the Lanczos results dis- 
play three different values of the local spin (see Fig.|SJ: A 
large and positive site, which corresponds to the strongly 
polarized "triplet" dimer, and two small sites (one pos- 
itive and one negative), building up the other "singlet" 
dimers. The local value of the magnetization as a func- 
tion of the spin-phonon coupling is reported in Fig. |UJ 
The spin texture has a rectangular shape, contains 12 
sites, and it is twelve-fold degenerate. 

As for the other cases, the Lanczos results for infinite 



FIG. 8: Spin density profile for 1/3-plateau. Full (empty) 
circles indicate sites with magnetization along (opposite to) 
the external field and the size of the circles is proportional 
to the spin amplitude. The bond lengths for a = 10 are also 
shown. 




FIG. 9: Local magnetization for the three different sites for 
the 1/3-plateau as a function of the spin-phonon coupling a. 



spin-phonon coupling coincide with the hard-core boson 
ones. As in the case of the 1/4-plateau, the spin tex- 
ture configuration shows diagonal "triplet" stripes, and 
this configuration is stabilized by the repulsion between 
strongly polarized dimers i^&SLSS* 2 ^ 



D. 1/2-plateau 

For completeness, we finally consider the calculations 
for the 1/2-plateau, by taking the 16-site cluster with 
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FIG. 11: Local magnetization for the two different sites for 
the 1/2-plateau as a function of the spin-phonon coupling a. 
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FIG. 10: Spin density profile for 1/2-plateau. Full circles 
indicate sites with magnetization along the external field and 
the size of the circles is proportional to the spin amplitude. 
The bond lengths for a — 10 are also shown. 

S% ot = 4. The typical spin configuration is shown in 
Fig. here, we have only two different sites, both with 
positive magnetization, corresponding to two different 
dimers: A "triplet" , strongly polarized, and a "singlet" , 
with a small magnetization. Hence, we obtain a square 
unit cell that contains only four sites. The evolution 
of the two local polarizations as a function of the spin- 
phonon coupling is reported in Fig. 1111 In this case, only 
the symmetry which interchange the "triplet" with the 
"singlet" is broken, whereas the translational symmetry 
is preserved. Therefore, the ground state is only two-fold 
degenerate, corresponding to the two possible choices of 
the "triplet" position in the unit cell. 



IV. COMPARISON WITH NMR 
EXPERIMENTS AT 1/8-PLATEAU 

As stated previously, at present, due to the very high 
magnetic field and the low temperature needed to stabi- 
lize the plateaus, it is only possible to perform accurate 
NMR experiments for the 1/8-plateau, which can give 
useful insight into the local spin textured Therefore, in 
the following, we make a more detailed analysis of our 
numerical results for this magnetization plateau. 

The Cu NMR spectra was measured at 35 mK in a 
field of 27.6 T, corresponding to the 1/8-plateau. The 
overall shape of the spectra can be well reproduced as- 
suming at least 11 distinct sites [see Fig. ll2f e)]. If we con- 
sider only the on-site dominant hypcrfine coupling, the 
hypcrfinc field is written as H n = A c g c (S z ), where (S z ) 



is the time-averaged local magnetization. The coupling 
constant A c and g c - value are determined as A c = —23.8 
T/fJ.B and g c = 2.28 by electron spin resonance (ESR) 
and NMR measurements Thus, the 11 sites observed 
in NMR indicate the existence of 1 1 distinct spin sites in 
the 1/8-plateau. Two large positive (S z ) sites (two neg- 
ative hyperfinc field sites) and one large negative (S z ) 
site (one positive hypcrfine field site) can be read off 
from Fig. 1121' e') . In addition to them, several sites spread 
around zero. 

In this plateau, two types of unit cells (square and 
rhomboid cells) have been proposed theoretically and so 
we calculated the spin textures in both unit cells as dis- 
cussed in Section IIII Al In both cases, there are two 
large positive spin sites and one large negative spin site, 
which hardly depend on the shape of the unit cell. How- 
ever, the distribution of the expectation values of the 
spin component (S z ) on the other sites is different in 
the two cases. In the square unit cell, the values (S z ) 
are concentrated around zero. On the other hand, in 
the rhomboid one, the expectation values spread further 
away from zero. These facts indicate that the results on 
rhomboid cell are qualitatively consistent with the results 
of NMR. Hypcrfine fields assuming the on-site hyperfine 
coupling A = A c = —23.8 T//ib are shown in Fig. I12f a) 
(Fig. I12f c)) for square (rhomboid) unit cell, in the case 
of a = 10. Following Ref. l32l the hyperfine fields includ- 
ing the effects of transfer hyperfine couplings B and C 
(see Fig.^) with the square (rhomboid) unit cell are also 
shown in Fig. ll2f b') fFig. H2f d)). Here the parameters B 
and C are chosen to reproduce the two large positive sites 
reported in the experiments. These results support the 
realization of the rhomboid cell. However, it is difficult 
to estimate the realistic values of the transfer hypcrfinc 
couplings and therefore the agreement of the rhomboid 
cell with the experiment is only indicative. 

It seems that the superstructure in the rhomboid unit 
cell is qualitatively consistent with the state observed 
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FIG. 12: Histogram of the hyperfine field (H n ). (a) The 
spin distribution in the square unit cell assuming only the 
on-site hyperfine coupling A — A c — —23.8 T/{ib- (b) The 
spin distribution in the square unit cell assuming the transfer 
hyperfine coupling B and C. On-site hyperfine coupling is 
A — A c — AB — C. (c) The spin distribution in the rhomboid 
unit cell assuming only the on-site hyperfine coupling A = 
A c — —23.8 T/fiB- (d) The spin distribution in the rhomboid 
unit cell assuming the transfer hyperfine coupling B and C. 
(e) The hyperfine field observed in NMR measurements p— 



by NMR. There are still quantitative differences though. 
For example, the numbers of different spin sites are not 
consistent: in the rhomboid cell, there are 8 different 
spin sites (see Fig.[3J) and in the experiments at least 11 
sites exist. On the other hand, there are only 6 different 
spin sites in the square unit cell. In the following we will 
discuss the possible effects of the inter-layer coupling and 
of the Dzyaloshinsky-Moriya interaction on the local spin 
structure and we emphasize the possible role of inter- 
layer coupling to achieve a more quantitative agreement. 



A. inter-layer coupling 

The three-dimensional (3D) structure of SrCu2 (1303)2 
consists of CUBO3 layers, intercalated by magnetically 
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FIG. 13: The 3D orthogonal dimer model for SrCu 2 (B0 3 )2 



inert Sr layers. The magnetic ions Cu 2+ form a three- 
dimensional lattice structure shown in Fig. 1131 and we 
denote the inter-layer coupling by J". However, we ex- 
pect that antiferromagnetic coupling J" is much smaller 
than the intra layer interactions J and J' because of the 
presence of the Sr layers. It is worth noting that, for 
small J" /J, the product of the dimer singlet, Eq. J5J, 
is the ground state also for the 3D orthogonal dimer 
model4iii£ In addition, in this limit, the magnitude of 
the triplet excitations and their dispersion relation do 
not depend on the inter-layer interaction J" ml These 
facts further support the hypothesis that the magnetic 
properties of SrCu2 (1303)2 can be very well described 
by the 2D orthogonal dimer model. However, the inter- 
layer interaction J" may affect the magnetic properties 
at high temperatures or under external magnetic fields. 
For instance, the inclusion of J" is necessary to correctly 
reproduce the behavior of the magnetic susceptibility at 
high temperatures ii£ 

Therefore, the effects of the inter-layer coupling might 
affect the superstructures at the plateaus. Although it is 
not presently possible to perform an exact calculation by 
Lanczos diagonalization which includes a full 3D lattice, 
we try to combine the previous exact results with the 
perturbation theory based on the hard-core bosons to 
extract some conclusion on the possible 3D spin texture 
and on the possibility to obtain more than 8 sites inside 
the unit cell. Therefore, let us begin by considering the 
hard-core boson picture, and, following Rcfs. I2f3l2 : nl28l 
we calculate the inter-layer interactions between triplets 
in two neighboring layers (which we denote by A- and 
B-plane, see Fig. 1141 by using the perturbation theory 
in the limits J'/ J -C 1 and J" / J <C 1. Starting from a 
state where one of the singlets is promoted to a triplet 
on each layer, up to fifth order in J'/ J, the triplets are 
completely localized and the excitation energy is given 
by: 

E = 2A + W k , (6) 
where A is the spin gap energy for one triplet and W k 



B -plane 




(3) 
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FIG. 14: Inter-layer interactions between the triplet excita- 
tions. The nearest-neighbor interaction Wi, and the next- 
nearest-neighbor interactions Wi and W'% are shown. 



the interaction between the two triplets at distance k, see 
Fig. 1141 Up to third order, the spin gap energy is written 
as A = J[l - (J'/J) 2 - l/2( J'/J) 3 ] and the interaction 
energies are given by: 
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The interaction energies vanish for k > 2. Note that, sim- 
ilarly to what happens in the third-neighbor intra-laycr 
interaction, the interactions for k = 2 strongly depend 
on the relative position of the triplets ,2§i2L2£i 

In the following, we consider the case of the rhomboid 
unit cell, which seems to be a better condition to repre- 
sent the real compound. We remind that, in this case, 
we have 8 different sites, and therefore in each unit cell 
there are two equivalent sites. From the hard-core bo- 
son calculation, we can assume that the triplets on the 
B-layer are located as far as possible from the triplets on 
the A-layer. In this assumption, the stacking pattern for 
the 1 /8-plateau is given in Fig. El When the inter-layer 
coupling is considered, the two triplets of the cell repel 
each other and the two equivalent dimcrs in the cell of 
the A-layer face two different dimers in the cell of the 
B-layer. Therefore, because of the inter-layer couplings, 
the two equivalent sites in a given plane can split into 
two different local spins: Assuming the stacking pattern 
of Fig. El we expect 14 different spin states per unit cell. 
In general, if we consider long-range interactions between 
two neighboring planes, we can obtain up to 16 different 
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FIG. 15: A stacking pattern for two neighboring layers. Black 
circles are in the A-plane whereas white circles are in the B- 
plane. The numbers indicate the different spins within the 
unit cell (see Fig. [3J ; the primed numbers are used to distin- 
guish the spins modified by the inter-layer interaction. 



spins per unit cell. However, the splitting for the sites 
indicated by 1 and 2 in Fig.El ma y be negligible because 
the difference originates from long-range interactions. 

Therefore, from this analysis, it come out that the in- 
clusion of the inter-layer coupling can give rise to 14 — 16 
different spin sites, which is in closer agreement with the 
NMR results, that indicate at least 11 different sites per 
unit cell. Note that in the square unit cell the number of 
different spin sites remains 6 even including the effect of 
inter-layer coupling J" . 



B. Dzyaloshinsky-Moriya interaction 

Finally, we consider the effect of the Dzyaloshinsky- 
Moriya interaction on the spin texture. Indeed, this term 
is relevant for SrCu2 (603)2: ESR and inelastic neutron 
scattering experiments show an anisotropic behavior of 
the spin gap, which depends on the direction of the ex- 
ternal fieldi 3 ^ 3 - and such a behavior can be explained by 
considering the Dzyaloshinsky-Moriya interaction^ 3 ^ 

If, following Ref. I4H , we assume in first approximation 
that the CUBO3 layer as the mirror plane, there is a 
Dzyaloshinsky-Moriya term only on the J' bonds and its 
component is perpendicular to the plane. In that case, 
the Hamiltonian for the Dzyaloshinsky-Moriya coupling 
is 



Wdm — D z 



Oj Oj J. 



(10) 



(n.n.n.) 



Note that Dzyaloshinsky-Moriya term 7^dm is odd under 
the exchange i <-> j and therefore, we have to fix the 
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FIG. 16: Local magnetization for the eight different sites 
in the rhomboid unit cell for the 1/8-plateau as a func- 
tion of the spin-phonon coupling for different values of the 
Dzyaloshinsky-Moriya interaction D z . 



direction from i to j for a pair i, j, as shown in Fig. ^ 

Strictly speaking, in the real compound, there is a 
slight buckling of the CUBO3 planed and so other higher 
components of Dzyaloshinsky-Moriya interaction can ap- 
pear. However, the magnitude of this buckling is very 
small and therefore the higher terms might be neglected 
with respect to the one given by Eq. ^| In practice, 
we have included the Dzyaloshinsky-Moriya interaction, 
Eq. I|10[) . in the original Hamiltonian and we have 
performed our self-consistent Lanczos diagonalization of 
the 2D cluster. The results for the rhomboid cluster 
with D z /J = 0.02 and 0.04 are shown in Fig. 1161 and 
compared with the ones for D z = 0. The inclusion 
of the Dzyaloshinsky-Moriya interaction makes the uni- 
form state more stable for small spin-phonon coupling, 
and the inhomogencous state appears for a spin-phonon 
coupling larger that the one we obtain without this in- 
teraction. The reason for this shift is due to the fact 
that the Dzyaloshinsky-Moriya term favors the hopping 
of the triplet. On the other hand, for strong enough spin- 
phonon couplings, i.e., a > 7, the values of the local spins 
are only slightly modified by the Dzyaloshinsky-Moriya 
interaction, and the number of different spins is always 
equal to 8. 

Therefore, we can infer that the Dzyaloshinsky-Moriya 
interaction is not a fundamental ingredient in determin- 
ing the spin texture at the magnetization plateau, and, 
in a first approximation, can be neglected. 



V. CONCLUSION 

In this paper, we have considered the effect of adi- 
abatic phonons on the 2D orthogonal dimcr model by 



using a self-consistent Lanczos diagonalization of small 
clusters. This study is directly related to the proper- 
ties of SrCu2 (1303)2, a new spin gapped material, under 
magnetic field. From an experimental point of view, at 
present, the only available informations about the mag- 
netic texture are about the 1/8-plateau, which is stabi- 
lized for magnetic fields around 27.6 Tcsla. Our theoreti- 
cal results for this plateau indicate that the spin-phonon 
coupling is able to stabilize two possible candidates for 
the local spin texture: i) a square unit cell of 16 sites 
with 6 different values of the spins and ii) a rhomboid 
unit cell of 16 sites with 8 different values of the spins. 
Although the energy difference between the two configu- 
rations is very small, the rhomboid cell contains a larger 
number of different spins and therefore it describes bet- 
ter the experimental finding. A closer agreement with 
the NMR results can be achieved by the inclusion of 
the inter-layer coupling, that can split equivalent sites, 
whereas the Dzyaloshinsky-Moriya interaction does not 
seem to play an important role in determining the actual 
values of the local spins. From our calculations, it comes 
out that the superstructure at the 1/8-plateau is quali- 
tatively consistent with the NMR results, but it is still 
open question to make a more quantitative comparison. 
The main limitation of our Lanczos calculations is the 
presence of finite size effects. As we already emphasized 
previously, the finite size effects might be small because 
the decay of the Friedel-like oscillation of the local spins 
is quite fast. However, in a small cluster, the uniform 
state is easily stabilized, and, therefore, the spin-phonon 
coupling required might be bigger than the realistic value 
(in the thermodynamic limit). Second, as mentioned in 
Section llVl the inclusion of the transferred hyperfine cou- 
plings to neighboring spins is probably needed to reach a 
better agreement. Indeed we showed that the inclusion 
of the transferred hyperfine coupling to the nearest- and 
next-nearest-neighbors improves the agreement between 
experiments and theoretical results. Moreover, it might 
be necessary to include also transferred hyperfine cou- 
plings from sites on neighboring layers. Assuming the 
stacking pattern in Fig. 1151 such hyperfine couplings can 
induce different hyperfine fields at two equivalent sites of 
a given plane. Thus these hyperfine couplings enhance 
the effects of the inter-layer couplings which we discussed 
in Section llVAl However, it is difficult to obtain reliable 
estimates of the transferred hyperfine couplings, and this 
discussion can be only indicative. 

Although at present, there arc no experimental results 
on the higher plateaus, we studied also the effects of the 
spin-phonon coupling for 1/4-, 1/3-, and 1/2-plateaus, 
with the hope that in the near future it will be possible 
to have experimental insight into the local magnetiza- 
tion of these plateaus (or at least some of them). Our 
results indicate that stripe-like superstructures are real- 
ized at the 1/3- and 1/4-plateaus, whereas a structure 
with a square unit cell is stabilized at the 1/2-platcau. It 
is worth noting that all the results regarding the shape 
of the unit cells arc consistent with the hard-core boson 
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picture, although within this hard-core boson approxima- 
tion it is not possible to obtain more than two different 
sites. In this respect, our Lanczos results give valuable 
insight into the real magnetization pattern that appears 
in the fascinating regions of the magnetization plateaus. 
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